
S0 = 20000;
mu = 0.03;
sigma = 0.4000;
K = 20000;
T = 30 / 365;
N = 100;
M = 100000;

eur_call(S0, mu, sigma, K, T, N, M)
eur_put(S0, mu, sigma, K, T, N, M)

underlyingpath_dw(S0, mu, sigma, T, N)

asian_call_ctrl(S0, mu, sigma, K, T, N, M)
asian_put_ctrl(S0, mu, sigma, K, T, N, M)

S0 = 20000;
mu = 0.03;
sigma = 0.4000;
K = 20000;
T = 30 / 365;
N = 100;
M = 100000;
b_ratio = 0.01;

barrier_put(S0, mu, sigma, K, T, N, M, b_ratio)

ceshi(S0, K, sigma, mu, T, N, M)